Phase Diagram and Spectroscopy of FFLO states of two-dimensional d-wave 

superconductors 
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Experimental evidence suggests that the Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) state may be 
realized in the unconventional, heavy-fermion superconductor CeCoIns . We present a self-consistent 
calculation of the field versus temperature phase diagram and order parameter structures for the 
FFLO states of quasi-two-dimensional d-wave superconductors. We calculate the spatially nonuni- 
form order parameter, free energy density, and local density of states for magnetic fields parallel to 
the superconducting planes. We predict that the lower critical magnetic field transition between 
the spatially uniform and nonuniform FFLO state is second order. We discuss the signatures of the 
nonuniform FFLO state which should be observable in scanning tunneling microscopy measurements 
of the local density of states. 
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I. INTRODUCTION 

The Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) state is 
predicted for clean spin-singlet superconductors as a re- 
sult of the competition between pairing correlations, fa- 
voring anti-parallel spin alignment, and the Zeeman ef- 
fect, favoring parallel spin alignment along the fieldi^ 
The compromise is a spatially inhomogeneous state of 
"normal" and "superconducting" regions. The "nor- 
mal" regions are defined by a spectrum of spin-polarized 
quasiparticlcs. The high-field FFLO phase was origi- 
nally suggested for superconductors with fcrromagncti- 
cally aligned impuritics^S but it was soon realized that 
a FFLO state should develop in superconductors in an 
external field if the Zeeman coupling dominates the or- 
bital coupling. 3 ! 4 ! 5 

A number of layered organic superconductors have 
been suggested as candidates for FFLO phases ALSiSiiSiii 
However, recent interest has focussed on the heavy- 
fermion material CeCoIn5^ This material has quasi- 
two-dimensional (2D) metallic planes, and shows evi- 
dence of unconventional, d-wave, superconductivitjti 3 * 1 ^ 
with a transition temperature of T c = 2.3 K. The quasi- 
2D tetragonal crystal structure consists of alternating of 
Ccln3 and Coln2. Superconductivity is believed to de- 
velop on a Fermi surface sheet which is nearly cylin- 
drical, i.e., a sheet in which the conduction electrons 
move primarily in the 2D metallic planesii 5 - CeCoIns 
exhibits the necessary characteristics for the existence 
of the FFLO state, ft is in the clean limit as indi- 
cated by the small normal-state residual resistivity and 



the small Sommerfeld coefficient for the low-temperature 
specific heat in the superconducting stated The up- 
per critical field is paramagnetically limited, as indicated 
by the large ratio of orbital critical field to Pauli criti- 
cal fields and CeCoIns exhibits a critical point where 
the transition from the normal to superconducting state 
changes from second-order to first order below the tem- 
perature T = 0.7 or 1.3 .fiT, depending on the field 
orientationiMiiSiiLiSiiSiSS Finally, several different exper- 
iments indicate a second-order phase transition within 
the superconducting state, which is assumed to be the 
transition into the FFLO state. 20 ! 21 ! 22 ! 23 ! 24 ! 25 

For s-wave superconductors it has been predicted that 
the upper critical field of the FFLO state depends on 
the dimensionality of the Fermi surfaced Fermi liq- 
uid interactions^ orbital effects and the presence of 
vorticesj£i2& and impuritiesA 2 ^ 3 ^ The one-dimensional 
stripe pattern for the order parameter originally pro- 
posed by Fulde and Ferrell, and Larkin and Ovchinnikov 
has been generalized to include higher harmonics and 
more complicated 2D FFLO latticesi 26 ! 31 ! 32 Burkhardt 
and Rainer included Fermi liquid effects and examined 
self-consistcntly quasi-2D superconductors^ while Mat- 
suo et al. formulated their analysis for three-dimensional 
(3D) Fermi surfaces^ In 2D the normal to supercon- 
ducting transition at B C 2 is second order, unless Fermi 
liquid effects are large. In 3D the FFLO transition is 
first order at high temperatures, but second order at suf- 
ficiently low temperatures^ For both 2D and 3D cases, 
the transition at the lower critical field, B c \, from the 
uniform superconducting state to the non-uniform FFLO 
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state is second order and signaled by the appearance of 
a single domain wall. Generally both spin and orbital 
coupling to the field are present. This leads to coexis- 
tence of the FFLO with vortices; these states have been 
studied extensively and it is predicted that the resulting 
inhomogeneous state has a number of novel properties, 
including a change in the order of the transition as well as 
additional phase transitions within the nonuniform phase 
corresponding to transitions between states belonging to 
different Landau levelsi 6 i 28 . 34 i 35 i 36 i 37 . 38 i 39 i 40 

For d-wave superconductors the FFLO phases are 
modified by the intrinsic anisotropy of the order param- 
eter. In the class of quasi-2D superconductors the upper 
critical transition line, B C 2(T), has a kink at low tem- 
peratures corresponding to the discontinuous change in 
direction of the stripe modulation with respect to the 
nodes of the order parameter ^^Mi^ However, the re- 
orientation transition, and the kink, are absent in 3D^ 
or if the shape of the Fermi surface is changed^ Re- 
cent calculations of the spatial modulation of the or- 
der parameter in 2D near B C 2 predict that the energet- 
ically favored state at low temperature and high field 
forms a "square lattice" instead of stripe order ^24^ The 
interplay between Pauli paramagnetism and the vortex 
phase of d-wave superconductors has also been studied 
in Refs. JHH). 

Theoretical predictions for the properties of the FFLO 
phases depend on material properties, including the 
quasiparticle interactions. Accurate results for the ther- 
modynamic and transport properties of the FFLO state 
generally require self-consistent calculations of the inho- 
mogeneous order parameter structure, local excitation 
spectrum and quasiparticle self energy. This is often a 
formidable task as one must examine a great number 
of different, but energetically nearby, states. Most of 
the theoretical work mentioned above was limited to the 
vicinity of the upper critical field, or near the maximum 
critical transition temperature, T FFLO , of the FFLO state. 

In this paper we address the stability of the nonuniform 
FFLO phases of 2D d-wave superconductors over the field 
range from the lower critical field, B&, to the the upper 
critical field, B C 2- Our analysis is for applied fields par- 
allel to the superconducting layers. In this geometry the 
effect of the magnetic field on the superconducting con- 
densate enters only through the Zeeman coupling of the 
quasiparticle spin to the field. For simplicity we assume 
a cylindrical Fermi surface, which is also supported by de 
Haas- van Alphen measurements on CeCoInsi 1 ^ 

The outline of this paper is as follows: In Sec. ILTI we 
introduce the theoretical model and the relevant aspects 
of the quasiclassical formulation of the equations for in- 
homogeneous superconductors. The B-T phase diagram 
for two-dimensional d-wave superconductors is discussed 
in Sec. IIIII This is followed by Sec. IIVI where we report 
calculations of the local density of states in the FFLO 
phases. We examine the energetics of the stripe phases 
and address the question of how to observe the spatial 
modulation of the order parameter structures in Sec. [V] 



Conclusions are provided in Sec. IVII 

II. THEORETICAL MODEL 

The quasiclassical theory of superconductivity 4 ^ 8 ^ 4 ^ 
is well suited for the study of spatially inhomogeneous 
states of superconductors varying over distances large 
compared to the Fermi wavelength, e.g., the supercon- 
ducting coherence length, £o = u//(27rT c ) ^> kj 1 . We 
consider spin-singlet superconductors for which the or- 
der parameter factorizes, A(R, p) = A(R)y(p), into 
a spatially dependent complex amplitude, A(R), and 
a momentum-space basis function, y(p), for an even- 
parity, one-dimensional representation of the crystal 
point group, where p is the direction of the Fermi mo- 
mentum. The basis function is normalized, (y 2 (p)) F s = 
1, where (. . .) FS represents the angular average over the 
Fermi surface weighted by the relative angle-resolved 
normal-state density of states on the Fermi surface. 

In the following we consider the FFLO states of a clean 
d-wave superconductor with d x 2_ y i symmetry described 
by the one-dimensional basis function, 3^(p) = \/2 (p 2 ~ 
py) = V2 cos2</>p. In this case the magnitude of the zero 
temperature gap parameter is A /(27rT c ) = 0.241i^i 

We solve the Eilcnbcrgcr equation in a magnetic field 
B for the quasiclassical (Matsubara) Green's function 
5(R,p;e m ), 

[i£ m % - A - v, g] + iv f ■ V g = , (1) 

self-consistently with the order parameter, A(R). In ad- 
dition, the Green's function must satisfy the normaliza- 
tion condition 

f = -n 2 l. (2) 

The quasiclassical Green's function in Nambu space can 
be decomposed into scalar and spin vector components 

~_ ( 9 + S-cr (/ + f ■ er)«cr 2 \ 

v> 2 (.f + f'-<x) -g + g-cr* J ■ W 

Correspondingly, the order parameter matrix in Nambu 
space for a spin-singlet superconductor is 

The Zeeman coupling of the quasiparticle spin with mag- 
netic field is given by 

v = + { fjS ■ <7* J ' ^ 

where Ui are Pauli spin matrices and (i = {g/2)^B is the 
absolute value of the magnetic moment of a quasiparticle 
with negative charge e; /is = |e|/(2mc) is the Bohr mag- 
neton. Note that the ^-factor is a material parameter, as 
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FIG. 1: The phase diagram of the FFLO state of a two- 
dimensional d-wave superconductor. The Larkin-Ovchinnikov 
(LO) state is stabilized in the high-B, low-T region of the 
phase diagram. The solid lines are second order phase tran- 
sition lines that determine the upper critical field, B C 2, and 
separate the normal and FFLO states below Tfflo ~ 0.56 T c . 
Below T < 0.06 T c a first order transition (long-dashed line) 
occurs between order parameter modulations along (110) and 
(100) directions. At the lower critical field, B c \, a second 
order transition (circles-solid) occurs between the uniform 
and nonuniform (llO)-oriented LO phase. The unphysical 
transition line from the uniform state into the (lOO)-oriented 
nonuniform state is shown for comparison (short-dashed). 
The Chandrasekhar-Clogston phase transition line between 
the uniform superconducting and normal state would be of 
first order below T F flo (dot-dashed), but is unphysical, while 
it is of second order and physical above Tfflo (thick dot- 
dashed). The corresponding order parameter modulations 
for a field scan between points (a) through (d) are shown 
in Fig. 0. 



is the Fermi velocity, which incorporates the high-energy, 
short-wavelength renormalizations of the bare electron 
g-factor into the effective g-factor defining the Zccman 
coupling of the quasiparticle to the field. 

For a spatially uniform order parameter, A = A(p), 
the physical solution of the Eilenbcrgcr equation for the 
Green's function <?(p; e m ) is 
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where 



ie± = ie m =F e B 



D±(p) = V|A( P )| 2 -(ze±) 2 , 



b. 



(6) 



and b is the direction of the uniform magnetic field. 

The uniform solution is easily generalized to the Fulde- 
Fcrrell (FF) state in which the order parameter has the 



FIG. 2: Order parameter modulations at T = 0.15 T c for dif- 
ferent points in the LO state as indicated by (a)-(d) in Fig. 
The period of the order parameter decreases with increasing 
field B. At the lower magnetic field, B c \, a transition from the 
uniform to nonuniform superconducting state occurs at posi- 
tion (a), which is signaled by a single domain wall (kink). We 
define the superconducting coherence length £o = W//(27rT c ) 
and dimensionless magnetic field b = /j,B/(2nT c ), with the 
zero-temperature order parameter Ao/(27rT c ) = 0.241. 



spatially varying phase, A(R, p) = A(p) exp(iq • R). 
The solution in this case is obtained by replacing the 
frequency in Eq. © by ie± — > ie± — rj. where r\ = 
ivy • q. The Green's function is then obtained from 
the uniform solution by a uniform gauge transformation 
<7(R,p;e m ) = Ug(p;e m ) W where U = exp(i q • Rf 3 /2). 
The linearized form of these solutions are used to obtain 
the upper critical field, B C 2(T), for the FFLO states, as 
well as the order parameter and thermodynamic func- 
tions for the homogeneous phase below B c i(T). 

Obtaining the Green's function for a nonuniform order 
parameter is more challenging, and generally requires nu- 
merical solution of Eilcnbcrgcr's equation. Schopohi^iiS 
and Nagai^^ introduced an efficient and numerically 
stable method for solving Eilenberger's equation by 
transforming it to a Riccati-type equation. We use the 
Riccati method as described by Eschrig^ From the solu- 
tions to the Riccati equation we can construct the com- 
ponents of the Green's function. The diagonal Green's 
functions, g and g, determine the quasiparticle excita- 
tion spectrum, while the off-diagonal Green's functions, 
/ and f , determine the order parameter through the BCS 
gap equation. For pure spin-singlet d-wave pairing, 



c e m >0 J 



(8) 



/(R,p;£ m )+/'(R,p;e ro )*-2 



7rA(R)y(p) \ 

km | J 



The gap equation is solved self-consistently together with 
the Riccati equations. 

There are often many nearly degenerate order param- 
eter configurations. In order to examine the relative sta- 
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FIG. 3: Second order instability lines between normal and 
FFLO state as function of relative orientation <f>o of one- 
dimensional stripes, i.e., order parameter modulations. For 
temperatures above T* ~ 0.06 T c , the modulation along (110) 
has the highest instability field, while below this temperature 
the orientation (100) is energetically favored. The thin-dotted 
line is the unphysical critical field, or "supercooled transi- 
tion", if one assumes a second order transition between the 
normal and uniform (q = 0) superconducting state. 
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FIG. 4: Dimensionless wave vector Q — q£o of order pa- 
rameter modulation along the instability lines in Fig. 
Left panel: T-dependence of Q for different orientations, with 
(j>0 = 0° <-> (100) and 4> = 45° <-> (110). Right panel: Angu- 
lar dependence of Q (dots) and upper critical field B C 2 (solid 
line) at T = 0. 



III. 



PHASE DIAGRAM 



bility of different order parameter solutions, we also cal- 
culate the free energy obtained from a quasiclassical re- 
duction of the Luttinger-Ward functional 

F S {B,T)-F N {B,T) =/"dRA/(R), (9) 

where the difference functional is 

A/(R) =U' dAr WS TrS (M»)- 

(10) 

N F is the density of states per spin at the Fermi level of 
the normal state. Note that g\ is an auxiliary propagator 
obtained from the solution to the transport equation with 
the physical order parameter scaled by the dimensionless 
coupling parameter, < A < 1, 

[i£m% - AA - v, gx] + iv f ■ V g x = . (11) 

The zero-temperature value of the free energy of a uni- 
form 2D (i-wave superconductor differs from the one of 
an s-wavc superconductor not only by the angular Fermi 
surface averages, but also by an extra term due to the 
nodal regions of the order parameter, 

AF(T = 0) = -N F ( (|A( P )|2)fs - (12) 

-iV F (e B Vet-|A(p)|2 Q(e 2 B - |A(p)| 2 )) FS . 

Consequently, the Pauli limiting field of an isotropic s- 
wave superconductor is fj,B P = 1.25 T c . While for a d- 
wave superconductor /i£> P = 1.19 T c , due to the nodal 
regions on the Fermi surface and in particular due to the 
second bracket in Eq. I|12|) . 



Various aspects of the phase diagram for the FFLO 
states of a 2D d-wave superconductor, particularly the 
region near B C 2, are discussed in the literature. We ex- 
tend those studies down to the lower critical field, B cl . 
We consider one-dimensional order parameter modula- 
tions resulting in one-dimensional stripes. Neglecting 2D 
modulations could result in a modified phase diagram at 
ultra-high fields and ultra-low temperatures^ 

The form of the order parameter amplitude for the 
FFLO states is then 

A(x,p) = A(p)f(x), (13) 

where f(x + X) = f{x) is periodic with A = 2ir/q, and x is 
the projected distance along q; x — q-R/g. The function 
f(x) can be cither real or complex and is a generalization 
of the Fulde-Ferrell state with a complex order parameter 
solution, 

A FF (R,p) = A(p) e ^' R , (14) 

or the Larkin-Ovchinnikov state with a real solution, 

A LO (R,p) = A(p) sin(q-R). (15) 

In the vicinity of B C 2, the gap equation Eq. JHJl can be lin- 
earized, if one assumes that the transition into the normal 
state is of second order. Therefore the LO state, whose 
order parameter solution is a linear superposition of FF 
solutions with wave vectors ±q, has the same instability 
line B C 2 as the FF state. 

The wave vector q, or the period A, are implicit func- 
tions of temperature and magnetic field. In the s-wave 
case all directions arc equivalent and the modulation, q, 
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FIG. 5: Angle- averaged, spin-dependent local density of 
states (LDOS) for a single domain wall solution. (Top): Pro- 
file of order parameter at T/T c = 0.15 and b = 0.175. The 
panels below show LDOS at locations (a) and (b). (Center): 
LDOS at position (a), where order parameter is zero, i.e., 
at domain wall. The Andreev bound states are shifted from 
zero energy by +6 for spin-up and — b for spin-down electrons. 
(Bottom): LDOS at position (b), far from the domain wall, 
i.e., in bulk of d-wave superconductor. The Andreev states are 
bound to the region where the order parameter is suppressed 
and are very weak far from the domain wall. 
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FIG. 6: Angle-averaged LDOS for a periodic solution inside 
the LO state. (Top): Profile of order parameter at T/T c = 
0.15 and b — 0.200. The panels below show LDOS at locations 
(a) and (b). (Center): LDOS at position (a), where order 
parameter vanishes. The broadened Andreev bound states, 
quasiparticle resonances, are shifted from zero energy by +b 
for spin-up and —b for spin-down electrons. (Bottom): LDOS 
at position (b) at maximum of order parameter. The Andreev 
bound states decayed and broadened even more. 



can point along any direction in real space. Thus, the 
translational invariancc is broken spontaneously in the 
FFLO state. In a d-wave superconductor the intrinsic 
anisotropy of the order parameter in momentum space fa- 
vors specific directions for the broken translational sym- 
metry. When q points along the nodal direction of y(p), 
we refer to this FFLO phase as the (110) stripe orienta- 
tion, and when q is along the anti-nodal direction it is 
the (100) orientation. 

Finally, in order to calculate the phase diagram we 
need to compute self-consistently the order parameter 
and find the free energy as a function of the triplet of 
parameters, T, B, X. The physical FFLO state that is 
realized at a given temperature and field, (T,B), cor- 
responds to a period A that minimizes the free energy. 
The resulting phase diagram of FFLO states is shown in 
Fig. p. 

We find that the real (Larkin-Ovchinnikov) order pa- 
rameter is favored over the complex (Fuldc-Fcrrcll) order 
parameter everywhere in the nonuniform FFLO region of 
the phase diagram. The upper critical fields, £? C 2, for dif- 
ferent stripe orientations shown in Fig. ^ (solid and con- 
tinuing into dotted lines) cross at T* w 0.06 T c , resulting 
in a kink in the upper transition line, in agreement with 
previous worki22i2ii2iii£ On the other hand, the lower 
critical field, B&, for the stripe orientation along (110) 
(line with circles) is always lower than the critical field 
for the stripe orientation along (100) (short-dashed line). 
This indicates a first order transition between states with 



different orientation of stripes (long-dashed line). Except 
for a small region in the FFLO phase diagram at ultra- 
high fields and ultra-low temperatures, where the (100) 
phase becomes favorable, the spatial modulation along 
the (110) direction is more stable. 

The order parameter profiles for a series of magnetic 
fields in the LO phase are shown in Fig. [3 As in the s- 
wave case, the period of the FFLO state decreases with 
increasing magnetic field. The transition between the 
uniform superconducting and FFLO states is second or- 
der and is signaled by the appearance of a single domain 
wall. The lower critical field transition is analogous to 
the s-wave scenario^*^ which contradicts the claim of a 
first order transition by Yang and Sondhy4£ The upper 
critical field transition from the FFLO to normal state 
is also second order with the amplitude of the order pa- 
rameter vanishing continuously with increasing magnetic 
field. 

Figure shows the second order instability line from 
the normal state into the nonuniform FFLO state for dif- 
ferent orientations of the stripes with respect to the anti- 
nodal direction, denoted by the angle (f>o- This transition 
is obtained by linearizing Eqs. and JHJ with respect 
to the order parameter A(R, p) = Ao3^(p) exp(iq ■ R). 
One of two special directions of q for the order parame- 
ter modulation, i.e., (100) and (110), yields the highest 
instability fields for all temperatures. 

In the left panel of Fig. 01 the dimensionless instability 
wave vector Q = q £o of the nonuniform LO state is plot- 
ted as a function of temperature for the different transi- 
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FIG. 7: Spin-up/down DOS averaged over a single period 
of the order parameter for three different b fields in the LO 
state. A broad bound state at e/(2ttT c ) — b is seen at lower 
fields, but continues to broaden as b increases. Finally, the 
DOS becomes flat and "normal" -like. Here T/T c = 0.15 and 
the corresponding upper critical field is b C 2 = 0.25. 



tions shown in Fig.|3 The right panel of FigEJshows the 
zero-temperature variation of the critical field B C 2 and 
wave vector q as a function of the stripe orientation <j)Q. 
The evolution of the upper critical field and instability 
wave vector with stripe orientation and temperature is 
the basis for the phase diagram shown in Fig. ^ 



IV. QUASIPARTICLE DENSITY OF STATES 

Regions in which the order parameter changes sign, as 
in the LO phase, have an associated spectrum on low- 
energy quasiparticles, which in this case is also spin- 
polarized. The local quasiparticle density of states 
(LDOS) is easily calculated from the retarded Green's 
function <7 R (R, p;e), which is determined from Eq. (JTJ 
using the sclf-consistently determined order parameter 
and analytic continuation to real energies from the upper 
half plane: ie m — > s+i0 + . The spatial and angle-resolved 
densities of states (DOS) for spin-up and spin-down ex- 
citations along the magnetic field are defined by 



N ul) (R,p;e) 



1 



3m (g R ± g R • b) 



(16) 



In Figs. 15181 we show the evolution of the LDOS and spa- 
tially averaged DOS for spin-up and spin-down excita- 
tions for T/T c = 0.15 as a function of magnetic fields 
and position in space. The spin-down DOS are shown 
as solid lines, while the spin-up DOS are indicated by 
dashed lines. 

Fig-Elshows the anglc-avcragcd, spin-dependent LDOS 
for a single domain along the nodal (110) direction. This 
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FIG. 8: An example of a spatially averaged DOS for (100) 
modulated order parameter at T = 0.15 T c and b — 0.200. 
Note that the physically stable solution has (110) modula- 
tion. (Upper left): Profile of periodic order parameter with 
three marked positions. (Lower left): DOS at the domain wall 
(x = 0), the Andreev bound state is clearly seen at e = e B . 
(Lower right): DOS at distance x = A/8 from domain wall, 
the Andreev bound state is slightly suppressed. (Upper right): 
DOS at the maximum of order parameter, the Andreev bound 
state is completely suppressed. In addition to the broadened 
Andreev bound state at e = e B there are singularities due to 
multiple Andreev reflections. 



configuration is in many ways equivalent to the scenario 
of the order parameter suppression near a (110) interface 
in a clean d-wave superconductor, c.f. Barash et ali££ 
There is a strong Andreev bound state at the position 
of the domain wall, arising from the change of the sign 
of the order parameter along quasiparticle trajectories 
crossing the domain wall (middle panel). In a magnetic 
field the Andreev bound state is no longer pinned to the 
Fermi level, but is shifted to e = e B = fJ-B (—/j,B) for 
spin-up (spin-down) quasiparticles. As we move away 
from the domain wall the bound state decays and the 
bulk DOS is eventually restored (lower panel), albeit with 
weak resonances at e = ±e B . 

Andreev bound states and resonances are characteris- 
tic features of the FFLO state. They exist as a result of 
sign changes of the order parameter, and are visible in 
the DOS for the periodic (110) stripe phases as well, see 
Fig. The resonance broadens due to hybridization of 
the energy levels associated with adjacent domain walls. 
Consequently the sharp bound state for a single domain 
wall broadens into a band that fills the low energy re- 
gion below the maximum gap to the continuum spec- 
trum. The quasiparticle resonance is visible as a broad 
peak centered at e = ±e B , even in the spatially aver- 
aged DOS as shown in Fig. [7| Note also the evolution of 
the spectrum with increasing magnetic field (from bot- 
tom to top panel). As the field increases the resonance 
broadens and the DOS approaches the normal-state value 
N(e)/Nf = 1, as is expected for a second-order transition 
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FIG. 9: (Upper left): Profile of order parameter for a single 
domain wall along (110). (Lower left): Profile of free en- 
ergy density A/ for single domain wall. (Upper right): Pro- 
file of order parameter for periodic domain walls along (110). 
(Lower right): Profile of free energy density A/ for periodic 
domain walls. 
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FIG. 10: (Upper left): Profile of order parameter for a single 
domain wall along (100). (Lower left): Profile of free en- 
ergy density A/ for single domain wall. (Upper right): Pro- 
file of order parameter for periodic domain walls along (100). 
(Lower right): Profile of free energy density A/ for periodic 
domain walls. 



to the normal state. 

Evidence of the FFLO state could be obtained from 
scanning tunneling microscopy by direct observation of 
the Andreev resonances organized along one-dimensional 
stripes, similar to the spectroscopic observation of vortex 
arrays in type II superconductors^ Since the quasipar- 
ticle resonance spectrum and distribution in space de- 
pend on the magnetic field, these states should be eas- 
ily differentiated from resonances arising from scatter- 
ing by impurities or interfaces ,SfliSiiS2i For the same rea- 
sons, point-contact spectroscopy, similar to recent mea- 
surements performed on CeCohi^^^A but for magnetic 
fields in the FFLO phase, should be able to identify the 
characteristic signatures of the FFLO state. 

In Fig. [S] we show the DOS for stripes oriented 
along (100). In contrast to the (110) orientation, 
there is a much richer spectrum along the (100) direc- 
tion. Additional singularities are visible at e/(2irT c ) ~ 
0.02,0.12,0.28,0.38. These singularities originate from 
quasiparticles that propagate nearly parallel to the 
stripes. For these trajectories there is a broad dip in the 
order parameter amplitude that produces bound states 
due to multiple Andreev reflections from the "walls" of 
this "order parameter well" . These resonances gener- 
ally occur at finite energy, in contrast to the topologi- 
cal bound states which are pinned to zero-energy in zero 
field. The bound states are further shifted by the field, 
and the topological bound states - shifted by the Zeeman 
field - arising from trajectories crossing the stripes are 
also visible in the spectrum. 

The spectrum for (100) stripes is similar to the LO 
state in s-wave superconductors because quasiparticles 
propagate along trajectories with similar order parame- 
ter profiles for directions perpendicular and parallel to 
the stripes. However, the anisotropy of the d-wave gap 



function in momentum space leads to more singularities 
below the maximum gap compared to the isotropic s- 
wave case. 

Finally, we note that Maki, Won and 
collaborators^^*^ calculated the spin-averaged DOS 
for the (100) striped FFLO state near the upper critical 
field at T — 0. Their result for the spatially averaged 
DOS, expanded to second order in the order parameter, 
is qualitatively different from our sclf-consistently 
calculated DOS at T/T c = 0.15 and b = 0.2 (see Fig.®; 
sufficiently so that a direct comparison is not possible. 



V. FREE ENERGY DENSITY 

The bound state and resonance spectra associated with 
the FFLO phases shown in Figs. [S] and [S] for the quasi- 
particle density of states also lead to structure in the free 
energy densities. For the stripe orientation along (110) 
see Fig. |5J and for (100) stripes see Fig. ^3 The free 
energy density exhibits additional structure for stripes 
oriented along (100) compared to (110). At the domain 
walls the free energy density is that of the normal state, 
while it is minimum where the magnitude of the order 
parameter reaches its maximum. The overshooting be- 
havior of the free energy density within a few coherence 
lengths away from the domain wall, as seen for a single 
domain wall in the lower left panels of Figs. |5J and I1UI 
accounts for the energy gained compared to the uniform 
superconducting state. 

The free energy density plots for a single domain wall 
shown in the lower left panels of Figs. 151 and 1 101 illustrate 
the subtle differences in energy between order parame- 
ter modulations along (110) and (100) directions. Only 
a self-consistent calculation of the nonuniform order pa- 
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ramcter can resolve such detail. 



VI. CONCLUSIONS 

We presented self-consistently calculated phase dia- 
gram of the FFLO state of a two-dimensional d-wave su- 
perconductor within the quasiclassical theory. We found 
that the Larkin-Ovchinnikov state is favored over the 
Fulde-Ferrell state in two dimensions. The oscillations 
of the order parameter amplitude along the nodal (110) 
direction for the gap function cos 2<f) is stable over most of 
the FFLO phase diagram. However, for one-dimensional 
stripes the state with (100) orientation is stable in the 
ultra-high B region of the phase diagram and below 
T* w 0.06 T c . Both B C 2 and B c \ transitions are of sec- 
ond order. The order parameter vanishes continuously 
at the upper critical field B c2 . At the lower critical field 
B c i a single domain wall, similar to a Bloch domain wall, 
signals the onset of the nonuniform FFLO state. 

The local density of states calculations suggest that 
the Andreev resonance spectrum may be used to iden- 
tify experimentally the intrinsic structure of the FFLO 
phases. In particular, the topological Andreev bound 
state is shifted from zero energy by an applied mag- 
netic field. All the quasiparticle resonances live in re- 
gions where the order parameter is suppressed. Thus 
their spatial distribution depends on the wave length of 
the nonuniform modulations, which are controlled by the 
strength of the magnetic field. Besides the character- 



istic energy dependence of the resonances, their spatial 
periodicity in the FFLO state enables us to tell them 
apart from random impurities in scanning tunneling mi- 
croscopy measurements. 

Finally, we note that the phase diagram calculated 
for the model two-dimensional <i-wave superconductor 
does not resemble the experimental phase diagram of 
CcCoIns. Thus, the identification of the observed phase 
transition in CeCoIns with a FFLO phase of a c?-wave 
superconductor is cither incorrect, or inconclusive and 
a more detailed model for the FFLO phases in this sys- 
tem needs to be implemented, e.g., a model with includes 
impurity scattering effects, Fermi liquid effects (e.g., ex- 
change interactions), orbital magnetization effects as well 
as more realistic Fermi surface topology. All of these ma- 
terials effects arc known to result in changes to the order 
and the location of the phase transition lines. 
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